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We consider the statistics of time delay in a chaotic cavity having M open channels, 
in the absence of time-reversal invariance. In the random matrix theory approach, 
we compnte the average valne of polynomial fnnctions of the time delay matrix Q = 
-ihS^dSjdE, where S is the scattering matrix. Our results do not assume M to 
be large. In a companion paper, we develop a semiclassical approximation to S- 
matrix correlation functions, from which the statistics of Q can also be derived. 
Together, these papers contribute to establishing the conjectured equivalence between 
the random matrix and the semiclassical approaches. 
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I. INTRODUCTION 


Quantum scattering processes at energy E can be described by the scattering matrix 
S{E), which transforms incoming wavefunctions into outgoing wavefunctions. This matrix 
is necessarily unitary, in order to enforce conservation of probability and, consequently, 
conservation of charge. Another important operator is the Wigner-Smith time delay matrix 
a hermitian matrix related to the energy derivative of S' as Q = -ihS^dSjdE. Its 
eigenvalues are the delay times of the system, and its normalized trace is the Wigner time 
delay, Ty/ = These quantities contain information about the time a particle spends 

inside a scattering region. A thorough discussion can be found in the review-^. 

We consider a scattering region (‘cavity’) inside of which the classical dynamics is strongly 
chaotic, connected to the outside world by small, perfectly transparent, openings. This 
can be realized in experiments with microwave cavities^^-, quantum dots^“— and compound 
nuclei^. We assume M open channels, so that S and Q are M-dimensional. We also assume 
there is a well dehned classical decay rate F, such that the total probability of a particle to 
be found inside the cavity decays exponentially in time as The quantity td = l/F is 
called the classical ‘dwell time’. 

When the wavelength is much smaller than the cavity size, the S and Q matrices are 
strongly oscillating functions of the energy, and a statistical approach is advantageous. One 
such approach is based on random matrix theory (RMT). Its main hypothesis is that S 
behaves like a random unitary matrix, distributed in the unitary group U{M) according to 
some probability measure (in the presence of time-reversal invariance, S must also be sym¬ 
metric; we do not consider that situation in this work). For ideal openings, this distribution 
is the normalized Haar measure of the group. We denote the averages with respect to this 
distribution by the symbol (•). In particular, the average of the Wigner time delay is equal 
to the classical dwell timei^iM, (r^y) = td- 

This RMT approach has had much success in describing transport statisticsl^^— , such 
as conductance, shot-noise, etc. It can also be applied to time delay statistics, but usually 
not starting from the S matrix, but rather from the Hamiltonian of the system. This allows 
better control of the energy dependence and calculation of correlation functions, but requires 
mapping the problem to a nonlinear supersymmetric cr-modeU2i^— . 

On the other hand, Brouwer, Frahm and BeenakkeriS^ have found the joint probability 
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distribution for the eigenvalues of Q, let us denote them by This allowed the 

calculation of marginal distributions^, distribution of Wigner time delay (for M = 2,^ and 
in the limit M » 1,^) and the average value of moments^i^ 

1 M 

v(„ = ^Tr[(3"] = yTr, (1) 

A few other polynomial functions have also been computed^i^. A recent review, also 
considering extension to non-ideal openings and other symmetry classes, can be found inSi. 

We hereby advance the RMT approach to statistics of time delay, obtaining an explicit 
formula for the average value of general polynomial quantities of the kind 

( 2 ) 

for any hnite set of positive integers ni,n 2 ,... Our method starts from the result oS^ and is 
based on Schur function expansions and determinant evaluations. Importantly, our results 
are not perturbative in the number of channels, being valid at hnite values of M. 

In the next Section we briehy present our results. Section 3 contains some preliminary 
material, and in Section 4 we present our calculations. 


II. RESULTS 


The average value of the moments have been found for general number of channels 
but expressed as a sum with M terms. Our results imply the following simple general 
formula, which contains a more efficient sum, with only n terms: 


where 


IK, X -kY 


[xY = x(x + l)---(x + n- l), [x]n = x(a;-l)---(x-n + l), 


( 3 ) 

( 4 ) 


are the raising and falling factorials. 

The hrst four cumulants of the Wigner time delay have been computed using some non¬ 
linear differential equation for their generating functionJ^ This amounts to hnding the value 
of YY) for j up to 4. Our results imply the explicit general formula 

[MY 


US') = ^ E s 


Ai-n 


[M]Y 


( 5 ) 
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where the sum is over all partitions of n, the length of a partition A is denoted i{X) and 


£(A) £(A) 

[Mf = n [^ - * + [M]x = n [^ + ^ - 1 ]a. 

2=1 2=1 


( 6 ) 


are generalizations of the rising and falling factorials. The quantity dx is the dimension of 
the irreducible representation of the permutation group labeled by A, and it is given by 


e{x) 

dA = n! n 


(A.-z + £(A))! AA, 


£(A) 

n (Xi-Xj-i + j). 


(7) 


The above examples are derived from particular cases of our most general result, which 
is the average value of a general Schur function of Q: 

( 8 ) 

nl [M\x 

These functions are homogeneous symmetric polynomials in the eigenvalues of Q. Since any 
symmetric polynomial in these variables can be expressed as a linear combination of Schur 
functions, this can be seen as a complete solution to the problem of computing the average 
value of polynomials (or analytic functions, if we allow inhnite series) in Q, such as the 
quantities Afni,n 2 ,... dehned in (|2]). For instance, the hrst of these which are neither of the 
form (E]) nor of the form (E]) are 

2M2(M2 + 2) 

“ (M2-l)(M2-4)’ 

and 



IM,A= 4M2(M^.8M2-3) __ 

\ 2,2/ \ 3.1/ (M2-1)(M2-4)(M2-9)' 

III. PRELIMINARIES 
A. Partitions and permutations 


A weakly decreasing sequence of positive integers, A = (Ai, A 2 ,...) is called a partition of 
n, denoted by A 1 - n or by |A| = n, if K = n. Each of the integers is a part, and the total 
number of parts is the length i{X). 

Partitions of n label the conjugacy classes of the permutation group S^- the cycle type 
of a permutation tt is a partition whose parts are the lengths of the cycles of vr, and two 
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permutations 7r,cr have the same cycle type if and only if they are conjugated, i.e. if there 
exists r such that vr = rar”^. Let C\ denote the set of permutations with cycle type A, and 
\C\\ the number of elements in C\. 

For any hnite group, there are as many irreducible representations as there are conjugacy 
classes. Therefore, partitions of n also label the irreducible representations of The trace 
of permutation tt, in the representation labeled by A, is denoted as and called its 

character. The character of the identity, xxi^) = d\, is the dimension of the representation, 
for which there is the explicit formula ([7]). Characters of Sn are class functions, i.e. X\i'^) 
depends only on the cycle type of tt and we may write X\{lA if tt € C^. Characters satisfy 
orthogonality relations, 

—^ 71 ^ 

E Xt.ir)xxira) = —Xx{(T)Sf,,x- (11) 

T€5„ «a 


B. Symmetric functions 


Let X be a matrix of dimension N, with eigenvalues Xi, 1 <i < N. Power sum symmetric 
functions of matrix argument are dehned as 


i(X) N 

P^{X) = nPA.(Jf). =Tr[Jf"] = yif. 

2=1 2=1 


( 12 ) 


They are clearly symmetric functions of the eigenvalues. 

Another important family of symmetric functions are Schur functions, related to power 
sums by^ 


‘>a(V) = - E \CAxx{p.)p^{X), p,(X) = E X,.(A)s„(V)- 

/il-72 /il-n 

These functions can also be written as a ratio of determinants, 

det(a;^'’~'^’^'^) 




where 


A(X) 


N N 


A(X) = det(x^ ) = n n - ^0, 


2=1 Jl=2 + 1 


is the Vandermonde determinant. 

The value of the Schur function when all arguments are equal to 1 is 

nl 


(13) 


(14) 


(16) 


( 16 ) 
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where [A^]^ is the generalization of the rising factorial dehned in (| 6 ]). Noticing that in the 
formula for d\ there appears the Vandermonde for Xi = Xi- i, it is also possible to show that 

N-l 

A({A.-i})=Sj(l«)nj!. (17) 

i=i 

Let dx = dxydxj^. In view of the identity 

J dxdet(Mxj)) det(g,{xj)) = N\ det (/ dxfi{x)gj{x)y (18) 

easily proved using the Leibniz formula for the determinant, the representation flT^ of Schur 
functions is useful for performing multidimensional integrals involving these functions. 

IV. STATISTICS OF THE TIME DELAY MATRIX 
A. Average of Schur functions 

We wish to compute the average value of a Schur function of the time delay matrix, s^(Q). 
This will be done using the following result obtained in2^: the probability distribution of the 
matrix 7 = is 

P( 7 ) = iA( 7)|2 det ( 7 )^ 6 -*^"^^^. (19) 

where 

n 00 

Z= / |A(7)Pdet(7)^e-^"°^^d7 (20) 

Jo 

is a normalization constant. 

Let Ti, 1 < i < M, be the eigenvalues of Q and 7 ^ = l/xi be the eigenvalues of 7 . The 
normalization constant is computed using Eq. fllSp : 

X oo A/f I 

d 7 det( 7 f+*-^e-^^^^^)det( 7 f^) = -^j^—^det{{M + j + i-2)\). ( 21 ) 

Standard determinant manipulations yield 

1 M 

The quantity we are after is 

1 

{sx{Q)) = — d7|A(7)|2det(7)^e-^^^^^SA(7'^)- (23) 

Z Jo 
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Writing the Schur function as a determinant, as in Eq. ffT^ . and using the following identity 
for the Vandermonde, 


A(7-') 


det 7^-1 


(24) 


we arrive at 


{sxm = ( rd7det(7r"*-'e-^"-^0det(7;"^"^'""). (25) 

Z Jo 


Using Eq. ffTSj) again we have 

MQ )) = + i - 2)!), (26) 

Consider the determinant det((a:j + i)!). Suppose we factor out a term (xj + 1)! from 
each row. The remaining determinant has the following structure: its ij element is a monic 
polynomial in xj of degree i - 1. It is well known^S that it therefore must be equal to the 
Vandermonde A (a:). 

Applying the above argument to fl26|l . we get 

1 M 

(*aW)) = 1)!, (27) 

where we used A({M - Aj + z - 2}) = - i}) and the special value of the 

Vandermonde, Eq. on. Plugging in the values of sa(I^) and Z, we arrive at 

= ( 28 ) 

or, in terms of the generalized falling factorial dehned in (|6]), our claimed result, 

(sa(Q)) = (29) 

n\ [MJa 


B. Particular cases 

The relation between power sums and Schur functions, Eq. flTSH . allows the calculation 
of more familiar quantities, such as 

(-^n) = ^(Pn(Q)) = E X\{n){sx{Q))■ (30) 

The character Xxi^) is different from zero only if X= (n- k, 1^) (so-called hook partitions), 
and is equal to (-1)^ in this case. On the other hand, the dimension d\ becomes (”'^^) for 
hooks, and with this we arrive at our example ([3]). 
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The other example we mentioned in Section 2 was 


(r^) = .i)(0)) — E d^MQ))- (31) 

Finally, consider the general moments A4ni,n2,...- We may assume that ^ = (ni,?7,2,...) is a 
partition of some integer, |/x|. Then, we have M.ni,n 2 ,... - and 

(Afni,n2,...) = E X\ip){sxiQ))- (32) 

Using this expression, we recover our examples (E]) and fITII]) . 
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